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Abstract 

We analyse the one-loop fermionic contribution for the scalar effective potential in the tem- 
perature dependent Yukawa model. In order to regularize the model a mix between dimensional 
and analytic regularization procedures is used. We find a general expression for the fermionic 
contribution in arbitrary spacetime dimension. It is found that in D = 3 this contribution is 



finite. 
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1 Introduction 



Recently there has been much interest in the phase structure of theories involving scalar fields 
presenting spontaneous symmetry breaking. Many aplications have been done in the Weinberg- 
Salam model and in grand-unified theories. The temperature generally is the parameter whose 
variation induces the transition from the broken to the unbroken phase, at least for the most 
current systems that develop first or second order phase transitions. 

To describe a second order phase transition the variation of the mass with the temperature is 
the most important fact. On the other hand the dependence of the coupling constant with the 
temperature may induce a first order phase transition in the scalar sector, as suggested by two of 
the authors in a recent work [I]. 

We start from the Yukawa model and we analyse the contribution coming from the fermionic 
loops for the temperature dependent scalar effective potential. The ultraviolet divergences are 
dealed with the method of analytic regularization || . We recall that the basic idea of this technique 
is to replace the denominator of the propagator (p 2 — m 2 + ie) by (p 2 — m 2 + ie) 1+a where a is 
the regulating parameter initially taken to be large enough. Consequently in a open connected 
set of points in the complex plane a the Feynman amplitudes are analytic. Then it is possible 
to analytically continue the Feynman expressions to the whole complex plane. In the Laurent 
expansions of these expressions we can identify the counterterms as the polar terms in the analytic 
extensions at some points of the complex plane. 

2 



To deal with finite temperature field theory using the imaginary time formalism we will have to 
use dimendional regularization in the momenta and deal with the Matsubara sums using another 
method. The most popular method to deal with the Matsubara sum is an analytic extension away 
from the discrete complex energies down to the real axis, with the replacement of the energies 
sums by countour integrals @ . If we are interested in systems at high temperature the decoupling 
theorem Q] allow us to use the dimension reduction method (DR). This approach has been used by 
many authors [||]. The basic idea is that in the imaginary time formalism the free propagator has a 
form (u>n+ p 2 + m 2 )^ 1 . The Matsubara frequency act like a mass so in the high temperature regime 
the non-static (n ^ 0) modes decouple, and we have a three dimensional theory. In other words, 
the only modes whose contribution do not fall of exponentially at distances much greater than (3 
are the (n = 0) modes of the bosons. Integration over the fermionic modes and the non-zero modes 
of the bosons result in a three dimensional theory. Of course this effective model will describe 
the original model only for distances R >> (3. As was stressed by Landsman || the standard 
summation method || based on analytic continuation do not work in the dimensional reduction 
approach. Instead, we have first to compute momentum integrals using dimensional regularization 
and to deal with the Matsubara sums, a inhomogeneous zeta function analytic regularization has 
to be performed. 

Recently such technique has been used to study different models at finite temperature. Ford and 
Svaiter [[?[ and Malbouisson and Svaiter [l] || studied the A(^ 4 and the Efimov-Fradkin (truncated 
or not) model at finite temperature. The possibility of vanishing the temperature dependent 
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coupling constants in these models has been investigated. In the first paper, assuming a non- 
simply connected spatial section, the thermal and topological contributions to the renormalized 
mass and coupling constant in the (Ac/? 4 ) d=a model was obtained at the one-loop approximation. In 
the second one, the authors extend the discussion of the massive self-interacting Ac/? 4 model to an 
arbitrary D-dimensional spacetime with trivial topology of the spacelike sections. The main result 
is that the possibility of a first order phase transition drived by the temperature dependent coupling 
constant, in the region where the model is super- renormalizable arises. The discussion in the case 
of a scalar model with non-polynomial interaction Lagrange density (the Efimov-Fradkin model) 
has been done in the third work. For D > 2 it was proved that at least two coupling constants 
of the truncated model may vanish and become negative by effect of temperature changes, while 
in the non-truncated model all the coupling constants remain positive for any temperature. The 
method used in the above quoted papers could provide an almost natural way to investigate 
stability regimes in finite temperature QFT models. 

It has been often sugested that the thermal contributions to the renormalized coupling con- 
stants of quantum models may bring up non-trivial effects. For instance, Gross, Pisarski and 
Yaffe @ argue that in finite temperature (QCD)^ the effective coupling constant g(A) decreases 
as the temperature or density is raised. In fact, they show in a perturbative context that at the 
first non-trivial order (QCD)^ should be asymptotically free at high temperature or pressure. In 
this approximation it is expected that at high temperatures thermal excitations produce a plasma 
of quarks and gluons which screen all (color) eletric flux. Such a transition from a low temper- 
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ature confined phase to a high temperature color screening phase has been also investigated by 



Polyakov [|K]] and Susskind |l Lj and others in lattice gauge theories. Such results have important 
astrophysical applications in the study of neutron stars or primeval universe models. 

The main goal of this paper is to investigate the one-loop femionic contribution to the scalar 
effective potential at finite temperature assuming that bosons and fermions interact via a Yukawa 
coupling. The outline of the paper is the following: in section II we briefly review the formalism 
of the effective potential. In section III the fermionic contribution to the effective potential is 
obtained. In section IV the singularity structure of the one-loop fermionic contribution to the 
scalar effective potential is studied Conclusions are given in section V. In this paper we use 
h — ks — c — 1. 



2 The effective action and the effective potential at zero 
temperature. 

In this section we will briefly review the basic features of the effective potential associated with 
a real massive self-interacting scalar field at zero temperature. Although the formalism of this 
section may be found in standard texbooks, we recall here its main results for completeness. Let us 
consider a real massive scalar field (f(x) with the usual \(p A (x) self-interaction, defined in a static 
spacetime. Since the manifold is static, there is a global timelike Killing vector field orthogonal to 
the spacelike sections. Due to this fact, energy and thermal equilibrium have a precise meaning. 
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For the sake of simplicity, let us suppose that the manifold is flat. In the path integral approach, 
the basic object is the generating functional, 

Z[J] =< 0,out\0,in >= 

J V[(p] exp{i[S[(p] + J d 4 x,J(x) V (x)]} (1) 

where V[ip] is an apropriate integration measure and S[(p\ is the classical action associated with the 
scalar field. The quantity Z[J] gives the transition amplitude from the initial vacuum |0, in > to 
the final vacuum |0, out > in the presence of some source J(x), which is zero outside some interval 
[— T, T] and inside this interval is switched adiabatically on and off. Since we are interested in 
the connected part of the time ordered products of the fields, we take the connected generating 
functional as usual. This quantity is defined in terms of the vacuum persistent amplitude 

by 

e iw[j] = Z[J]) (2) 
and the connected n-point functions G^(xi,x 2 , are 

r^(r r r ) - 6nw[J] i m 

Expanding W [J] in a functional Taylor series, the n-order coefficient of this series will be the 
sum of all connected Feynman diagrams with n external legs, i.e. the connected Green's functions 
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defined by eq.(3). Then 

00 1 r 

W \A = £ I / dW-d 4 x n G^{x 1 ,x 2 ....x n )J{x 1 )J{x 2 )..J{x n ). (4) 

n=0 n - J 



The classical field (fo{x) is given by the normalized vacuum expectation value of the field 

SW < 0,out\ip(x)\0,in >,j 
iPo[X) ~ 5J(x)~ <0,out\0,in>j ' [ ) 

and the effective action T[(p ] is obtained by performing a functional Legendre transformation 

r> ] = W[J\ - J d 4 xJ{x)ip {x). (6) 
Using the functional chain rule and the definition of tp given by eq.(5) we have 

= -J(x). (T) 

Just as W[J] generates the connected Green's functions by means of a functional Taylor expan- 
sion, the effective action can be represented as a functional power series around the value ipo = 0, 
where the coeficients are just the proper n-point functions r^ n \xi,x 2 , ~iX n ) i.e., 

00 1 r 

r [^o] = £ ~1 / d i x 1 d i x 2 ...d A x n T {n) {x 1 ,x 2 , ..,x n ) (po(xi)<p (x 2 )..(po(x n ). (8) 

The coefficients of the above functional expansion are expressed in terms of the connected one- 
particle irreducible diagrams {IP I). Actually, T^ n \xi,x 2 , ...,x n ) is the sum of all 1PI Feynman 
diagrams with n external legs. Writing the effective action in powers of momentum (around the 
point where all external momenta vanish) we have 

r> ] = J d*x(-V{^) + l -{d^fz[^\ + ..). (9) 



The term V{ipo) is called the effective potential [112] .To express V(ipo) in terms of the 1PI 
Green's functions, we write T^ n \x 1 ,x 2 , -.,x n ) in momentum space, 



T^\x 1 ,x 2 , ..,x n ) = — ^ J d 4 hd 4 k 2 ..d 4 k n (27r)*6(h + k 2 + ..k n ) e^ + - k " x ^ n \x u x 2 , ..,x n ). 

(10) 

Assuming that the model is translationally invariant, i.e. is constant over the manifold, we 
have 

I> ]= /d 4 a;f:l(f(")(0,0,...)(^or + ...)- (H) 

J n =l H - V 7 

If we compare eq.(9) with eq.(ll) we obtain 

%o) = -E3f (n) (0,0,..)br, (12) 

n n - 

then is the sum of the all IP I diagrams earring zero external momenta. Assuming that 
the fields are in equilibrium with a thermal reservoir at temperature in the Euclidean time 
formalism, the effective potential V(/3, (po) can be identified with the free energy density and can 
be calculated by imposing periodic (antiperiodic) boundary conditions on the bosonic (fermionic) 
fields. 

3 The one-loop effective potential of the Yukawa model 
at zero and finite temperature. 
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Let us consider a system consisting of bosons and fermions fields interacting via a Yukawa 
coupling in thermal equilibrium with a reservoir at temperature f3~ l . They are defined on a four 
dimensional flat spacetime with trivial topology of the spacelike sections. In the zero temperature 
case the generating functional for the scalar and fermionic fields correlation functions is given by: 

Z[rj,f],J] = J T>ipV'ipV(pexjp{i[S['ijj,'ip,(p\+ J d A xtpr] + fjip + J<p]} (13) 

where ip(x), ip(x), rj{x) and rj(x) are elements of the Grassmann algebra and <£>(x) and J(x) are 
commuting variables. 

The perturbativelly renormalizable action has the form, 

S®, if>, if] = J d A x Q(c^ 6 ) 2 - \ml V 2 b + V(<p b ) + Mi P - M - goVb^b) (14) 

where V(<fb) = jrjpi and mo, M are respectivelly the boson and the fermion bare masses and Ao 
and go are the bare coupling constants. Of course ipt and ijjy, are bare bosonic and fermionic field. 



The most general divergent terms are of the type f!3 | 



_ T div = j d A x (isZJdM 2 - Um 2 cp 2 + 

u " ' (15) 



+SZ^iip(i flip - 5Mifji(j - gSZgi/ji/jip + jdXcp 4 + ±5a(p 3 + 8apj 
Although the action is renormalizable, the model is not multiplicatively renormalizable. To 

circumvect this difficulty and to allow the theory to become multiplicatively renormalizable we shall 

introduce at the tree level action all terms which we expect to be generated by the renormalization 
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procedure, i.e. 



S(ip, -0, v?) = J d 4 x ^d/j,ipd^ip — ^ ^n<£ n + P — M — g^i/^j + counterterms , (16) 



where A 2 = |m 2 , A 3 = |, and A 4 = ^. 

As usual, perturbation theory is generated by, 



^3 c 

Z(ij, 77, J) = exp{-i J d^g-^j + V(—))}Z (fj, v , J) (17) 



where 



Z (fj, V , J) = exp-i J d 4 xd 4 y(r](x)A F (x - y)rj(y) + \j{x)A{x - y)J{y)) (18) 
with A F (x — y) and A(x — y) being respectively the fermionic and bosonic propagator functions, 

A F (x-y) = (i jd x + M)A(x-y), (19) 

and 



i r p-ip(x-y) 

A{x-y) = —- d 4 p- — . (20) 

(27t) 4 J p 2 — m 2 + xe 



From the above formulas, following a porcedure entirely analogous to that described in the 
preceding section for the pure scalar case it is easy to get the fermionic contribution to the 
effective potential V(ipo) 



10 



V(ip ) J (fx = iln J° dipdtp exp{i J d 4 xip(i jb — M — gcp)ip} (21) 

After a Wick rotation to Euclidean space and using the rules for Grassmann integrals we get 
the contribution from the single fermionic loops to the scalar effective potential 

V(<po) J d 4 x = -lndet(i jd E — M — gip ) (22) 
Using a well known result 



logdet(M + g<p ) = trlog(M + g<p ), (23) 

we have, 



lndet(i jd E -M- g<p ) = trlog{% jd E ) - J2 — —{M + gip Q ) s tr{— ) s . (24) 

Using a Fourier representation for and taking into account that the contributions from odd 
values of s in the above sum vanish, it is possible to recast the fermionic contribution to the 



effective potential in the form |I3 



V(9«) = ^J W? — e gy^. (25) 

In the finite temperature case using the Matsubara formalism we have to perform the replace- 
ments co — > uj n = jf{n + 5) an d ^ / dq E = 4 J2n- Then the contribution from the single fermionic 
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loops to the effective potential is given by 



vi n m 2 V V I (-!)' (M + gyp) 2 ' 



Let us define the quantities, 



and 



a=(±f (27) 



^> < 29 > 



where // is a parameter with mass dimension introduced to deal with dimensionless quantities 
performing analytic extensions. First we use dimensional regularization going to a generic D- 
dimensional spacetime. Then eq.(26) becomes 



\/(0,/?)=^f;aT- / ( As )( 7 + #r E - — ^7— t. (30) 

5=1 n=-oo {{n+ i) 2 ) S 2 

where /(D, s) is given by: 
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Before going one some commments are in order. It is well known Jl5| that dimensional reg 



ularization techniques for massless fields can not led to definite results due to the presence of 



infrared divergences fll5| . Since we are regularizing only a d = D — 1 dimensional integral, this 
procedure is equivalent to inserting a mass into the d dimensional integral. In other words, the 
Matsubara frequencies play the role of "masses" in the integral provided that we exclude the limit 
f3 — ► oo which means that we must restrict ourselves to non-zero temperatures. Another point is 
that in order to evaluate the one-loop finite temperature diagrams the usual approach is to express 
the integrand as a countour integral ||. In this paper we use another technique still aplying the 
principle of the analytic extension. 

In the next section we will analyse the singularity structure of the inhomogeneous Riemann 
zeta function and other factors appearing in eq.(30) in order to identify the divergent terms in the 
fermionic contribution to the effective potential. We start by analytically regularizing the model. 

4 The singularity structure of the fermionic contribution 
to the effective potential. 



As we remarked before the fermionic contribution to the effective potential is ill defined due 
to the singularities of the gamma function that appears in f(D, s) and the singularities in the 
Matsubara sum. The Matsubara sum may be expressed in terms of the generalized inhomogeneous 
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Riemann zeta function, which can be analytically extended to a meromorphic function in the whole 
complex s plane. The polar terms must be removed in the renormalization procedure. In order to 
identify these poles let us first recall the definition of the inhomogeneous Riemann zeta function 



or Hurwitz zeta function ITB 



OO 1 

«* J) =£<^p (32) 

which is analytic for Re(z) > 1. 

After some manipulations it is possible to express the Matsubara sum in eq.(30) in terms of 
C(z,q) and write V(4>,/3) in the form, 

Vfa (3) = 2^ D £ a*-f{D, s)( 7 + g<P) 2s ((2s - d, -). (33) 

8=1 A 

To analytically extend the inhomogeneous Riemann zeta function, we go along the following 
steps: first using the Euler representation for the Gamma function we write it as 

., -t 

POO p — 

C(z,q) = —- dtf- 1 - -. (34) 



r(z) Jo " 1 - e~* 

Next, we split the integral from zero to infinity in two integrals, from zero to one and from 
one to infinity. The second one is an analytic function of z, the divergences being associated to 
the zero limit of the first integral. Then using a Bernoulli representation for the integrand it is 
possible to get the following expression to the analytic extension of ((z, |) 
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where g%{z) is given by 



OO p 2 

2-1 C 



*W = rwA ?=T" <36) 



and the B n (x) are the Bernoulli coefficients [16]. We remark that in the literature there is 



another formula for the analytic extension of the inhomogeneous Riemann zeta function; the 



Hermite formula [17j given by 



1 q^ 2 f°° —z y 1 

C(z,q) = 1 \-2 (q 2 + y 2 )~ sinfz arctan -) — dy. (37) 

sv lH ' 2q z z-1 Jo V ; V q J e 2n y-l y V ; 

Of course the analytic extension must be uniquely defined and these are only different repre- 
sentation of the same analytic extension. Substituting the analytic extension given by eq.(35) in 
the fermionic contribution to the effective potential V((p, (3) we get, 



V(h0) = li D ET=i a%-h(D, s)( 7 + g<l>) 2a n _j± 1 +a+1) 

( f°° rli i^s-D _ei_ , v^°o gw(|) 1 
yjl ai 1 e t_i -T 2^ n =0 n \ 2s-D+n 

where the regular function h(D, s) is given by, 



(38) 



s T(s) 
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Let us analyse the two cases D = 3 and D = 4 separately. For the case D = 3 we have 

V(<f>, (3) = /i 3 ET=i « § -^(3, s) ( 7 + #) 2 W — 



r(-|+ s +i) 

(40) 

( poo J./, -i.25-3 e5 I y^oo 1 \ 

I Jl ««- << e *_i ~r Z^„=o n l 2s -3+n ^ 

The fermionic contibution to the effective potential is finite. There is no ultraviolet divergences 
in D = 3. One would not normally expect this since the tadpole graph is ultraviolet divergent 
(s = 1). This situation is very similar to the calculation of the renormalized vacuum energy 
of scalar fields confined in boxes (Casimir energy) JTB| . Dolan and Nash used the zeta function 
analytic regularization method to obtain the Casimir energy of conformally coupled scalar field 



confined in odd and even dimensional spheres [19]. They obtained that for odd dimensional 



spheres (even space-time dimension) there is a pole in the point of interested, being necessary the 
introduction of a counterterm, while for even dimensional spheres (odd dimensional space-time) 
the result obtained is naturally finite. No renormalization is needed. For a careful study of this 



subject see ref. [20 



For the case D = 4 we have 

VfaP) = // 4 Er=i« 2 ~ s M4,s)(7 + #) 2s ^ 



r(s-i) 

(41) 

roo ji .2s-4 , B o(k ) I jfeCgj _|_ V^oo gnCg) 1 

Jl aL <* e t_i -T 2s -4 " r 4s-4 " r ^n=3 n! 2s-4+n 

Note that the factor T~ 1 (s — 1) just cancels the pole from the term n = 2 in the sum over n. 
The pole comming from the term n = in the sum must be canceled by the introduction of a 
suitable counterterm. All other terms s > 3 are finite. 
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5 Conclusion 



The aim of this paper is to discuss an alternative method to deal with the Matsubara sum in a 
finite temperature field theory with bosons and fermions in interaction. We use this method to 
calculate the one-loop fermionic contribution to the scalar effective potential assuming the Yukawa 
coupling between fermions and bosons. Note that we are using a BPHZ scheme with subtraction 



at zero momentum of the Feynman integrals. Matsumoto, Ojima and Umezawa [[H]] claims that 
tha Matsubara method seems to produce temperature dependent divergences which disapear only 
after a summation over the Matsubara sums. We showed that the countertems are temperature 
independents. 

A curious observation is in order. We note that eq.(38) does not contains singularities for any 
odd space time dimension D, due to the fact that the sum is over integer values of s, and the 
Bernouilli coefficients -B n (|) — for n odd. For even values of D the fermionic contribution to the 
effective potential (see eq.(38)) has only a divergence due to the term s = n = 0. 

It would be interesting to generalize the method if we consider that there is a non-zero fermion 
density |22[|. This can be done introduzing a chemical potential a. At finite temperature the 



chemical potential will change the Matsubara frequencies by u n — > u n + ia [23|. In this case we 
have to analytically extend the inhomogeneous Epstein zeta function ((z, q) for complex q. This 
subject is under invetigation. 
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